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SUMMARY 


This  report  lays  the  groundwork  for  analyzing  the  feasibility  of 
producing  electromagnetic  signals  from  modulated  electron  beam  antennas 
in  the  ionosphere  and  detecting  the  signals  on  the  ground.  It  concen¬ 
trates  on  frequencies  between  tens  of  hertz  to  a  megahertz  because 
higher  frequencies  are  easily  radiated  from  conventional  antennas  on 
spacecraft.  This  study  is  geared  toward  helping  the  military  to  de¬ 
velop  low-frequency  communication  systems  that  use  beam  antennas. 

It  is  well  known  that  natural  electric  currents  in  the  ionosphere 
and  magnetosphere  produce  electromagnetic  radiation  of  frequencies  less 
than  1  MHz,  and  that  radiation  is  regularly  observed  at  the  ground. 
Certain  man-made  electron  beams  ejected  into  the  ionosphere  have  gen¬ 
erated  radiation  at  frequencies  above  1  MHz — this  radiation  has  been 
detected  on  the  ground.  As  yet,  however,  no  experiment  has  attempted 
to  use  modulated  electron  beams  in  the  ionosphere  to  produce  radiation 
of  the  modulation  frequency  or  its  harmonics,  and  to  detect  that  radia¬ 
tion  on  the  ground. 

This  detailed  and  extensive  analysis  necessarily  incorporates 
several  idealizations  that  tend  to  overstate  the  efficiency  of  the 
radiation.  Nonetheless,  this  report  suggests  that  a  carefully  con¬ 
trolled  experiment  may  be  able  to  produce  long-wave  radiation  that  can 
be  detected  on  the  ground. 

This  report  is  essentially  a  handbook  that  gives  predicted  char¬ 
acteristics  of  the  radiation  for  four  assumed  experimental  configura¬ 
tions.  The  first  is  hypothetical;  it  is  optimized  for  the  generation 
and  detection  of  long-wavelength  radiation  using  present-day  or  easily 
attainable  technology.  The  other  three  are  already  planned  and  funded 
(their  primary  emphasis,  however,  is  not  on  wave  generation  and  detec¬ 
tion  measurements) . 

Further  analysis  is  needed  to  design  the  optimal,  long-wavelength 
experiment  in  detail.  In  addition,  more  theoretical  work  is  suggested 
for  an  assessment  of  the  impact  of  relaxing  certain  of  the  idealiza¬ 
tions  made  in  the  beam/propagation  models  used  in  the  present  report. 
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This  work  represents  the  Pacific-Sierra  Research  Corporation  (PSR) 
final  report  under  contract  DNA001-84-C-0288.  The  research  was  spon¬ 
sored  by  the  Rome  Air  Development  Center  (RADC)  under  Air  Force  Office 
of  Scientific  Research  task  number  2305 J2,  and  the  contract  was  admin¬ 
istered  by  the  Defense  Nuclear  Agency  (DNA) .  The  task  performed  under 
this  contract  is  a  parameter  study  of  electromagnetic  emissions  from 
modulated  electron  beams  in  the  ionosphere.  It  will  be  used  to  define 
an  experiment  that  will  evaluate  the  feasibility  of  long-’avelength 
communication  using  modulated  electron  beam  antennas  in  the  ionosphere. 
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I.  INTRODUCTION 


This  Pacific-Sierra  Research  Corporation  (PSR)  report  is  a  param¬ 
eter  study  of  the  characteristics  of  electromagnetic  radiation  that  may 
be  emitted  by  modulated  nonrelativistic  electron  beams  in  the  mid¬ 
latitude  ionosphere.  Two  intentions  motivate  this  work.  The  first  is 
to  design  an  optimal  experiment  for  probing  the  physics  of  communica¬ 
tion  at  very  low  frequency  (VLF)  and  extremely  low  frequency  (ELF)  by 
means  of  an  electron  beam  as  a  virtual  antenna.  The  second  is  to 
examine  the  possibilities  for  communications  research  provided  by 
planned,  funded  electron  beam  experiments. 

Of  the  four  experiments  considered,  the  optimal  experiment  will 
employ  an  electron  gun  that  we  call  the  "future  gun."  Its  parameters, 
at  the  leading  edge  of  technical  feasibility,  will  provide  the  maximum 
possible  current  at  a  given  voltage.  The  future  gun  is  intended  for 
deployment  on  a  rocket  dedicated  to  communications  research. 

The  other  experiments  analyzed  here  are  not  dedicated  to  communi¬ 
cations,  but  they  will  utilize  electron  beams  that  are  injected  into 
the  ionosphere  and  that  are  or  may  be  modulated.  The  first  experiment 
will  be  flown  aboard  a  rocket,  BERT-I,  to  be  launched  in  April  1985. 
Devoted  mainly  to  spacecraft  charging,  BERT-I  will  use  a  relatively 
low-power  electron  beam.  Current  plans  are  for  an  unmodulated  beam, 
but  the  hardware  could  be  modified  to  provide  square  wave  modulation. 
The  second  utilizes  a  Vehicle  Charging  and  Potential  (VCAP)  experiment 
gun,  which  will  be  placed  on  Spacelab  2,  to  be  launched  in  mid-1985. 
Modulated  electron  beams  will  be  ejected  and  electromagnetic  fields 
will  be  measured  in  the  near-field  region.  The  third  experiment  is  a 
joint  Japanese  and  U.S.  venture  called  SEPAC  (Space  Experiment  with 
Particle  Accelerators),  which  will  be  placed  aboard  three  shuttle 
flights  in  August  1986.  The  gun  on  SEPAC  is  the  most  powerful  of  the 
three.  Its  beams  will  be  modulated,  and  near-field  measurements  are 
planned  for  it. 

The  mathematical  model  employed  for  all  the  calculations  presented 
here  is  from  Harker  and  Banks. ^  That  model  calculates  radiation 


emission  from  a  semi-inf inice  electron  beam  having  a  finite  correla¬ 
tion  length.  The  model,  which  is  reviewed  in  Sec.  II  of  this  report, 
is  coupled  with  simple  estimates  of  the  antenna  gain  and  atmospheric 
attenuation  to  provide  an  overall,  if  highly  simplified,  model  of  the 
radiation  received  on  the  ground. 

Section  III  encompasses  a  somewhat  detailed  analysis  of  the  future 
gun.  Because  the  interpretations  of  the  results  for  the  three  other 
guns  follow  the  patterns  established  for  the  future  gun,  those  analyses 
are  given  in  brief  and  only  a  few  figures  for  those  cases  are  presented 
in  this  report.  Sections  IV  through  VI  are  analytical  results  for  the 
planned  and  funded  guns — BERT-I  is  treated  in  Sec.  IV,  VCAP  in  Sec.  V, 
and  SEP AC  in  Sec.  VI. 


II.  THEORY 


MODEL 

The  Qarker  and  Banks'*'  model  assumes  chat  an  infinite  pulse  train 
of  electrons  is  emitted  from  the  electron  source  at  z  *  0  and  spirals 
around  a  magnetic  field  line  with  parallel  and  perpendicular  velocities 
V||  and  vj_.  The  beam  has  zero  diameter,  has  always  been  on  and  is  not 
modified  by  the  surrounding  plasma.  However,  the  magnitude  of  the 
current  falls  off  exponentially  in  z  with  an  e-folding  length  of  £. 

That  falloff  models  a  decreasing  ability  of  the  beam  to  radiate  with 
increasing  z. 

The  beam  passes  through  a  cold  plasma  in  an  ambient  magnetic  field. 
The  plasma  consists  of  electrons  and  one  ion  species  [which  we  assume 
to  be  (01^)+  in  this  parameter  study]  and  is  represented  by  the  relative 


dielectric  tensor 


and  e  ■  in  /m.  is  the  electron-to-ion  mass  ratio,  ft  ■  oj/oj  is  the  ratio 
e  i  ce 

of  the  (angular)  frequency  of  the  electromagnetic  wave  to  the  electron 

2  2  2  2 

cyclotron  frequency,  and  5  ■  (w  +  ui  . )/oi  is  the  normalized  sum  of 

pe  pi  ce 

the  squares  of  the  electron  and  ion  plasma  frequencies. 

Maxwell's  equations  are  solved  by  Fourier  analysis.  The  result¬ 
ing  formula  for  power  radiated  takes  the  following  form  in  the  case 
of  large  l: 


4*0*“.  s1”2 


3”— 00  i-1  p-1 


watts 


In  the  above,  c  is  the  velocity  of  light;  p  labels  the  harmonic,  so  that 

ui  “  pw„,  where  w  is  the  angular  modulation  frequency  of  the  beam; 
m  m 

i  labels  the  root  of  the  dispersion  relation;  and  s  labels  the  radia¬ 
tion  mode.  The  Cerenkov  mode  corresponds  to  s  -  0,  and  the  regular 

and  anomalous  cyclotron  modes  correspond  to  s  ■  1  and  s  *  -1,  respec- 

2  2 

tively.  The  quantities  y^  and  are  the  two  roots  or  branches  of  the 

dispersion  relation.  Properties  of  those  branches  are  summarized  in 

Sec.  4.3  of  The  Near-Earth  and,  Interplanetary  Plasma ,  an  excellent  book 
2 

by  Al'pert.  In  Eq.  (3),  1^  Is  the  beam  current  in  amperes,  RQ  is 

the  impedance  of  free  space  Rq  ■  377  ft,  b  is  the  length  of  the  electron 

pulsey  and  d  is  the  spacing  between  the  front  of  the  successive  pulses. 

2 

Finally,  the  quantity  is  the  square  of 


V1  T  2  2l1/: 

=1  -  -  ^  +  *  *3>  -  ViJ 

+  Ja(»)|  [(pi  -  -  El)  -  e'] 


where  J  and  J'  are  the  Bessel  function  of  order  s  and  its  derivative 
s  s 

and  n  is  the  parallel  index  of  refraction  n  •  U.  cos  0,  where  0  is  the 


angle  between  the  direction  of  the  ambient  magnetic  field  and  the  wave 
normal  (k  vector).  The  parallel  index  of  refraction  at  resonance 
satisfies  the  relationship 


and  the  quantity  u  is  the  gyroradias  times  the  perpendicular  component 
of  the  k  vector: 


The  solutions  are  further  limited  by  the  following  requirements: 

2  2 

(1)  y  must  be  real,  (2)  n  S  U  »  and  (3)  the  power  must  be  real  [i.e., 

the  product  of  any  two  of  the  radicands  in  Eq.  (4)  must  be  greater  than 

or  equal  to  zero].  Those  three  conditions  define  allowed  and  forbidden 
regions  for  emission  of  electromagnetic  radiation  from  a  modulated 
electron  beam.  The  boundaries  of  those  regions  are  complicated  func¬ 
tions  of  the  plasma  density  and  composition,  and  of  the  strength  of  the 
ambient  magnetic  field.  For  combinations  of  the  beam  energy  and  modu¬ 
lation  frequency  that  violate  these  requirements,  the  radiated  power 
is  set  to  zero  in  the  following  sections. 

The  e-folding  or  beam  correlation  length  l  in  the  above  is  neither 
known  nor  is  it  modeled.  In  this  report,  we  have  divided  out  the 
length  l  and  calculated  radiation  emitted  per  unit  beam  length.  A 

method  of  determining  the  minimum  beam  length  needed  to  receive  a 

signal  on  the  ground  under  the  most  favorable  circumstances  of  emis¬ 
sion  and  propagation  is  given  in  the  following  subsection. 

ROUGH  ESTIMATES  OF  PROPAGATION 

The  electromagnetic  radiation  from  a  modulated  electron  beam  will 
be  detectable  at  the  ground  if  its  power  at  the  receiver  is  higher  than 
that  of  the  background  noise.  Using  this  as  a  criterion,  one  can  esti¬ 
mate  the  minimum  beam  correlation  length  necessary  for  a  signal  to  noise 


ratio  greater  than  one.  The  following  analysis  is  based  on  the  linear 
theory  given  above. 

Implicit  in  Harker  and  Banks's'*"  theory  is  that  the  beam  antenna 
has  a  gain  G  that  depends  on  the  beam  correlation  length  l.  In  the 
case  of  a  beam  radiating  a  power  P  ,  if  the  receiver  is  a  distance  R 
from  the  electron  beam,  the  power/unit  area  P  that  reaches  the  receiver 
is  approximately 

P  “  “S  GP  ,  (7) 

47rtT  r 

where  A  is  a  transmittance  factor  that  depends  on  the  attenuation  of 
the  radiation  as  it  passes  through  the  ionosphere.  The  values  of  A 
reported  here  are  for  transmission  of  the  electromagnetic  radiation 
vertically  downward  and  ignore  possible  spreading,  refraction,  or  re¬ 
flection  of  the  electromagnetic  radiation.  The  ambient  received  noise 
in  a  particular  frequency  band  is  a  measured  quantity  often  expressed 
in  the  form  PqWq,  where  p^  is  measured  and  is  power/unit  area  per  hertz 
and  WQ  is  the  receiver  bandwidth  in  hertz.  The  criterion,  then,  for 
the  desired  detectable  signal  is: 

P  >  P0WQ  •  (8) 

In  Eq.  (7),  R  is  known  and  Pf  is  given  by  Harker  and  Banks^  in 
the  form 

Pr  “  Pr*  ,  (9) 

where  pr  is  the  radiated  power  per  meter  of  beam  correlation  length. 

The  transmittance  factor  for  radiation  through  the  entire  ionosphere 

3 

to  the  ground  surface  has  been  calculated  by  Booker,  Crain,  and  Field 
as  a  function  of  frequency  for  ambient  day  and  night  model  ionospheres. 

3 

We  are  justified  using  the  results  of  Booker,  Crain,  and  Field  for 
altitudes  greater  than  180  km,  because  most  of  the  attenuation  occurs 
below  180  km  altitude. 


The  beam  antenna  gain  G  is 


max,  radiated 


jl  r 

4ir  J  ■ 


?  ..  ,  <5Si 

radiated 


4^  V*  A(cos  x) 


or 


G  *“  A(cos  x)  ’ 


(10) 


where  x  la  the  ray  angle ,  pra<jiated  is  the  radiated  power,  and 

P  ...  is  its  maximum  value.  Equation  (10)  is  approximately 

max,  radiated 

correct  if  the  functions  in  the  integrand  of  the  power  expression 
[Eq.  (15)  in  Harker  and  Banks1  and  Eq.  (4)  above]  are  slowly  varying 
over  the  interval  where  the  term  [Eq.  (20)  from  Harker  and  Banks1] 

W(k||)  -  [(-/)  +  («  -  k!|V)1  -  suce)2J  ,  (11) 

is  appreciable.  Here  V|j  is  the  component  of  the  electron  velocity 
along  the  magnetic  field,  u  is  the  radiation  frequency,  ky  is  the  com¬ 
ponent  of  the  wave  vector  along  the  magnetic  field,  s  is  the  mode  number, 
and  u>ce  is  the  electron  cyclotron  frequency.  The  half  value  width  of 
the  function  W(k||)  is 


(12) 


Using  the  following  for  the  parallel  index  of  refraction: 


one  can  write 


A(cos  X) 


(13) 


where  0  is  the  angle  of  the  k  vector  with  respect  to  the  magnetic  field, 
and  Eq.  (12)  has  been  used.  The  functions  d6/dn  and  dx/d@  are  compli¬ 
cated  functions  of  the  ray  angle  X>  wave  normal  angle  9,  radiation  fre¬ 
quency  uj,  and  parallel  index  of  refraction  n.  The  function  d0/dn  is 
the  inverse  of  Eq.  (34)  found  in  Harker  and  Banks^;  dx/d0  is  defined 
in  Eq.  (38)  in  Harker  and  Banks. ^ 

Equation  (8)  can  now  be  written  as 


(14a) 


If  the  gain  predicted  by  Eqs.  (4)  and  (7)  is  unreasonably  high 

3 

(greater  than  300  or  10  ,  for  example) ,  experience  suggests  that  an 
upper  limit  be  artificially  imposed.  Equation  (7)  or  (14a)  then 
becomes 


4ttR2/pq\  1 
0’  A  \?JG0 


(14b) 
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Equation  (14a)  or  (14b)  gives  the  absolute  minimum  beam  correla¬ 
tion  length  2^  for  which  a  signal  may  be  detectable  at  a  distance  R 
away  from  the  electron  beam.  The  value  Zg  gives  a  lower  bound  to  the 
correlation  length.  Because  of  the  numerous  idealizations  in  the  beam 
and  propagation  physics  (all  of  which  ignore  effects  that  degrade  the 
signal),  lengths  greater  than  may  be  needed.  However,  beams  shorter 
than  lQ  will  not  suffice.  Therefore,  if  igis  short,  the  beam  probably 
can  radiate  sufficient  energy  for  detection  before  it  disintegrates, 
but  we  have  no  assurance  of  that.  For  combinations  of  parameters  that 
cause  Zg  to  be  very  long,  the  experiment  may  not  be  feasible  at  all. 

PRESENTATION  OF  CALCULATIONAL  RESULTS 

The  important  findings  are  displayed  in  graphical  form  in  the 
following  sections.  The  graphs  show  the  power  radiated  per  unit  length 
of  the  electron  beam,  plotted  against  the  kinetic  energy  of  the  elec¬ 
tron  beam  and  its  modulation  frequency.  Many  of  the  graphs  take  two 
forms.  The  first  form  is  a  three-dimensional  plot  showing  the  five 
highest  orders  of  magnitude  for  the  radiated  power,  which  is  useful 
for  easily  identifying  and  scanning  the  region  where  the  peak  power  is 
generated.  Figure  1  is  an  example  of  that  kind  of  plot.  The  second 
form  is  a  contour  plot  of  the  power,  which  is  useful  for  reading  numer¬ 
ical  values  of  the  power  and  for  identifying  null  or  forbidden  regions 
where  no  power  at  all  is  radiated.  Figure  2  is  an  example  of  the  con¬ 
tour  plots.  The  contours  in  all  the  figures  are  labeled  by  the  common 
logarithm  of  the  radiated  power,  in  watts  per  meter.  Note  that  the 
contour  labels  in  the  region  below  0.1  kHz  in  Fig.  2  are  preceded  by 
a  minus  sign,  whereas  above  10  kHz  many  contours  are  positive. 

In  some  regions  the  radiated  power  falls  quickly  to  zero  as  the 
modulation  frequency  or  the  beam  energy  changes.  These  regions  appear 
on  the  contour  plots  as  closely  spaced  parallel  lines  that  are  not 
labeled  because  of  lack  of  space  near  the  contours.  Such  a  region 
appears  on  Fig.  2  in  the  rectangle  between  10  and  1000  kHz  and  between 
60  and  100  keV.  On  some  of  the  figures,  small  unlabeled  diamonds 
appear  in  regions  otherwise  devoid  of  radiation.  These  indicate  the 
existence  of  narrow  regions  of  allowed  radiation  that  were  not  well 


mapped  fay  the  coarse  grid  calculated  for  the  contour  plot.  Such  regions 
are  not  of  interest  in  beam  experiments,  as  will  be  shown  in  Sec.  III. 

Each  figure  has  additional  labels  that  characterize  the  electron 

beam  and  the  ionospheric  plasma.  The  pitch  angle  of  the  beam  is  a  in 

degrees,  the  beam  power  is  P__.w  in  watts,  P  labels  the  radiation  har- 

monic  plotted,  s  labels  the  mode,  is  the  angular  electron  plasma 

frequency  in  radians  per  second,  ft  is  the  ratio  of  the  ion  plasma 

pi 

frequency  to  the  electron  plasma  frequency,  and  I_  in  amperes,  which 

o 

appears  on  some  plots,  is  the  beam  current. 

The  dispersion  relation  that  determines  the  magnitude,  direction, 
and  polarization  of  the  electromagnetic  radiation  has,  in  general,  two 
solutions  at  a  given  frequency.  Because  the  propagation  character¬ 
istics  of  the  two  solutions  differ,  each  is  treated  separately  in  most 
of  the  graphical  presentations  in  this-  report.  Each  curve  is  labeled 
"root  1"  or  "root  2"  according  to  whether  it  corresponds  to  one  or  the 
other  of  the  solutions.  Figures  1  and  2  are  root  1  solutions  and 
Figs.  3  and  4  are  root  2  solutions.  It  is  a  general  characteristic  of 
the  root  1  solutions  that  the  highest  powers  are  radiated  above  the 
lower  hybrid  frequency  and  below  the  electron  gyrofrequency  (from 
7.5  kHz  to  1.3  MHz  in  Figs.  1  and  2).  The  three-dimensional  and  con¬ 
tour  plots  look  like  a  sloping  mesa  or  plateau  in  that  region,  and 
for  that  reason  the  region  is  called  the  "central  plateau"  in  the 
sections  that  follow.  Figures  3  and  4  show  that  the  root  2  solutions 
tend  to  have  more  structure  than  the  root  1  solutions  and  fill  up 
regions  forbidden  to  root  1,  although  they  have  power  levels  generally 
much  lower  than  the  root  1  levels.  The  specific  trends  will  be  analyzed 
in  more  detail  in  the  following  sections. 
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III.  FUTURE  GUN 


CHOICE  OF  PARAMETERS 

The  parameter  studies  for  the  future  gun  were  based  on  the  observa- 
tion  that  the  most  powerful  nonrelativistic  electron  gun  flown  on  a 
rocket  to  date  delivered  30  kW  the  entire  duration  of  the  flight  (6  min) 
The  highest  current  achieved  so  far  is  about  100  A.  For  a  given  volt¬ 
age,  the  current  from  a  cathode  is  space-charge  limited,  a  relation  that 
is  often  expressed  in  terms  of  the  perveance  p  of  the  gun: 


P 


g 


IMAX 
v3/2  • 


(15) 


where  I  is  the  maximum  current  that  can  be  extracted  at  voltage  V. 
max 

Higher  currents  can  be  achieved,  however,  by  increasing  the  number  of 
cathodes,  changing  the  cathode  geometry,  or  by  using  several  guns  in 
parallel. 

It  is  not  clear  what  the  lower  limit  of  usable  gun  voltages  are, 
but  a  Japanese  sounding  rocket,  K-9M-41,  successfully  injected  thermal 
(0.1  eV)  electron  beams  into  the  ionosphere.^  Tsutsui,  Matsumoto,  and 
Kimura  did  not  give  the  beam  currents.  The  lower  limit  on  the  beam 
energy  is  limited  by  the  power  desired  in  the  beam.  To  obtain,  for 
example,  1  kW  in  the  beam,  a  10  V  beam  would  require  100  A  current. 

But  typical  electron  gun  perveances  are  on  the  order  of  a  few  times 
10  (A/V  ) ,  which  implies  that  10J  to  10°  cathodes  are  needed  for 
a  10  V  gun  to  deliver  1  kW  into  the  electron  beam.  Because  this  is 
an  unreasonable  number  of  cathodes,  and  although  we  have  no  reason  to 
suspect  (except  for  intuition,  based  on  terrestrial  communications 
experience)  that  1  kW  beam  energy  is  essential  to  communications,  we 
have  chosen  10  eV  as  the  lower  limit  to  the  beam  energy  in  our  analysis. 

The  upper  limit  to  achievable  voltages  is  in  the  relativistic 
range.  Our  analysis  is  inherently  nonrelativistic,  however,  so  its 
usefulness  is  questionable  much  beyond  50  keV.  Nonetheless,  we  chose 
10^  eV  as  the  top  end  of  the  parameter  study — not  because  we  deemed 


Che  results  at  such  high  energies  acceptable,  but  to  make  trends  stand 
out  more  clearly’ on  the  graphs.  As  the  pitch  angle  of  electron  injec¬ 
tion  and  the  electron  cyclotron  frequency  are  varied,  features  of  the 
radiated  power  shift  along  the  beam  energy  axis.  Thus,  power  peaks  that 
appear  in  the  nonrelativistic  theory  at  high  beam  energy  move  to  lower 
beam  energies  as  emission  parameters  change. 

The  analysis  of  the  future  gun  was  done  assuming  a  rocket  flight 
launched  at  White  Sands,  New  Mexico.  The  radiation  was  calculated  for 
electron  beam  emission  into  the  ionosphere  at  altitudes  of  100,  200, 

300,  and  400  km  for  day  and  night  flights.  Table  1  summarizes  the 
ionospheric  parameters  we  assumed  in  the  calculations.  We  chose  pitch 
angles  of  0.1,  45,  and  89.9  deg. 

One  set  of  calculations  was  performed  assuming  a  constant  beam 
power  of  30  kW  at  all  beam  voltages.  The  other  assumed  a  constant 
current  of  100  A  at  all  beam  voltages.  In  most  of  the  figures,  we 
show  only  the  power  that  goes  into  radiation  at  the  fundamental  fre¬ 
quency.  This  is  most  of  the  power  that  goes  into  the  radiation  for  a 
square  wave  modulated  beam  at  VLF.  It  Is  perhaps  1/10  to  1/50  of  the 
total  power  that  goes  into  the  radiation  at  lower  frequencies  for  a 
square  wave  beam  (the  rest  goes  into  higher  harmonics  in  the  allowed 
radiation  region).  For  a  sinusoidally  modulated  beam,  the  fundamental 
represents  all  the  radiated  power. 

Calculations  show  that  for  much  of  the  parameter  space  the  theo¬ 
retical  index  of  refraction  at  resonance  is  very  high  (in  some  regions 

6  9 

the  resonant  index  of  refraction  is  greater  than  10  or  10  )  for  radia¬ 
tion  into  the  regular,  anomalous,  and  higher-order  cyclotron  modes. 

In  those  regions  the  resonant  indices  are  much  larger  than  those  for 
radiation  into  the  Cerenkov  mode;  therefore,  total  internal  reflection 
of  the  signal  will  be  more  likely  for  cyclotron  modes  than  for  Cerenkov 
modes.  Radiation  in  the  Cerenkov  mode  may  therefore  dominate  the  radi¬ 
ation  that  actually  reaches  the  ground.  For  that  reason,  the  higher- 
order  modes  are  not  considered  except  in  a  sample  calculation  at 
200  km,  which  is  discussed  later. 

This  calculation  of  radiated  power  ignores  all  warm  plasma  and 
collisional  effects.  Those  effects  will  reduce  the  predicted  radiation 
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at  low  altitudes  (100  km)  and  will  reduce  the  predicted  resonant  values 
of  the  refractive  index  for  cyclotron  emission, 

MAGNITUDE  OF  RADIATED  POWER 

Figures  1  through  32  show  the  radiated  power  into  the  Cerenkov 
mode  at  the  fundamental  frequency  for  the  four  altitudes  considered 
here,  for  day  and  night  flights,  and  for  constant  beam  power.  In  all 
those  figures,  the  pitch  angle  has  been  set  to  45  deg.  Figures  1 
through  4  are  for  a  day  flight  at  100  km,  and  Figs.  5  through  8  are 
for  night  at  100  km.  Figures  9  through  12  are  for  a  day  flight  at 
200  km,  and  Figs.  13  through  16  are  for  night  at  that  altitude.  Simi¬ 
larly,  Figs.  17  through  20,  and  21  through  24  are  for  300  km,  day  and 
night,  respectively,  and  Figs.  25  through  28,  and  29  through  32  are  for 
400  km,  day  and  night,  respectively. 

Several  important  conclusions  can  be  drawn  from  these  figures. 
There  are  regions  in  parameter  space  where  there  is  no  radiation  at 
all,  and  the  locations  of  those  regions  change  with  altitude  and  de¬ 
pend  on  the  time  of  day.  The  easiest  way  to  get  an  overview  of  those 
dependences  is  to  look  at  contours  of  the  total  radiation  summed  over 
both  roots  of  the  dispersion  relation.  Figures  33  and  34  show  contours 
of  total  radiated  power  at  100  and  400  km,  respectively,  during  the 
day.  Figure  33  shows  that  at  100  km  during  the  day  there  are  three 
regions  in  which  there  is  no  radiation:  above  approximately  4  MHz  and 
less  than  300  keV;  between  10  kHz  and  2  MHz  and  above  100  keV;  and  a 
triangular  region  less  than  7  kHz  and  less  than  500  eV.  At  400  km, 
the  first  region  increases  in  size  and  joins  the  second,  so  there  is 
hardly  any  radiation  for  beam  energies  above  1.5  MHz,  and  none  from 
10  kHz  to  1  MHz  for  beam  energies  above  3  keV.  The  third  region  almost 
disappears  at  400  km,  though  a  remnant  of  it  exists  near  8  kHz  below 
100  eV. 

The  same  effect  occurs  at  night  (as  Figs.  35  and  36  show)  although 
the  limits  are  pushed  to  higher  beam  energies.  It  is  significant  that 
at  night  the  forbidden  region  at  the  low  frequencies  is  much  larger 
below  1  kHz  than  it  is  during  the  day. 

In  general,  the  structures  that  are  displayed  in  the  figures  move 
to  the  lower  beam  energies  as  the  altitude  increases.  (There  are  also 
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shifts  in  frequency,  which  are  somewhat  less  significant.)  The  values 
for  radiated  power  remain  more  or  less  constant  on  the  structures  as 
they  move,  although  the  radiated  power  at  a  given  frequency  and  beam 
energy  may  change  drastically. 

The  radiated  power  goes  mainly  into  root  1  on  the  central  plateau 

between  1  MHz  and  1  to  200  kHz,  depending  on  beam  energy.  The  values 
2 

range  from  10  W/m,  near  the  high  energy  edge  of  the  plateau,  to 
9 

10  W/m  at  low  energies.  The  power  into  root  1  below  about  50  Hz  is 

-5  -12 

much  lower,  typically  10  to  10  W/m.  The  radiated  power  into 

root  2  is  much  lower  than  into  root  1  almost  everywhere — typical  values 
-4  -16 

range  from  10  to  10  W/m  below  100  kHz,  In  part  of  that  region, 

however,  power  goes  into  root  2  only.  Above  1.5  MHz  for  daytime  emis¬ 
sion,  root  2  dominates,  particularly  at  low  altitudes  and  for  beam 
energies  less  than  10  keV. 

VALUES  OF  INDEX  OF  REFRACTION  AT  RESONANCE 

Detailed  computer  printouts  show  that  the  index  of  refraction  at 
resonance  usually  decreases  as  the  beam  energy  increases  and  as  the 
angle  of  the  wave  normal  moves  away  from  90  deg.  At  100  km,  the  index 

4 

of  refraction  is  a  maximum  of  about  10  at  10  eV  on  the  low-frequency 
side  of  the  central  plateau  of  root  1,  and  decreases  to  300  at  the  high 
frequency  side.  At  the  high  beam  energy  end,  the  index  of  refraction 
is  about  500  on  the  low-frequency  side,  and  about  6  on  the  high- 
frequency  side.  The  results  at  400  km  for  daytime  emission  are  simi¬ 
lar,  but  the  variation  at  the  high  energy  end  is  from  3000  at  the  low- 
frequency  side  to  300  at  the  high.  Figures  37  and  38  show  an  inter¬ 
mediate  case — altitude  200  km.  The  nighttime  values  are  similar. 

Figure  38  shows  that  the  peak  values  for  root  2,  around  2  MHz  for 
daytime  emission  at  low  altitudes,  vary  from  about  600  at  10  eV  to  1 
at  relativistic  energies.  The  printouts  also  show  that  values  for 
root  1  at  low  frequencies  (as  in  Fig.  37)  are  around  300  at  100  km  and 
1000  at  400  km  for  daytime  emission;  for  nighttime  emission  the  values 
range  from  10^  at  10  eV  to  40  at  relativistic  energies  and  at  low 
altitudes,  and  fall  around  600  at  all  energies  at  400  km. 
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WAVE -NORMAL  AND  RAY  ANGLES 


We  have  not  prepared  figures  showing  wave  normal  and  ray  angles. 
However,  detailed  computer  printouts  establish  that  the  angle  of  the 
k  vector  (wave-normal  angle)  and  the  ray  angle  vary  enormously  over 
the  parameter  space,  although  the  variation  with  altitude  and  time  of 
day  is  not  as  marked  as  that  for  the  index  of  refraction.  At  the  low- 
energy  end  of  the  root  1  central  plateau,  the  k  vector  ranges  from 
89.6  deg  at  the  low-frequency  side  to  34  deg  at  the  high-frequency 
edge.  The  corresponding  ray  angles  are  179.6  and  124  deg.  At  the 
high-energy  end,  the  k  vector  angle  ranges  from  89  deg  at  the  low- 
frequency  edge  to  42  deg  at  the  high-frequency  edge.  The  corresponding 
ray  angles  are  about  180  and  114  deg. 

The  wave-normal  angle  at  very  low  frequencies  is  close  to  90  deg 
at  relativistic  energies  with  a  ray  angle  also  close  to  90  deg.  At 
lower  energies,  the  wave-normal  angle  falls  to  between  20  and  30  deg 
at  low  altitudes  for  daytime  emission,  with  a  ray  angle  between  10  and 
30  deg.  For  nighttime  emission  at  low  altitude  and  daytime  emission 
at  high  altitude,  the  ray  angle  is  close  to  180  deg  for  root  1  and 
90  deg  for  root  2.  For  the  peak  in  root  2  at  low  altitude  and  2  MHz 
frequency,  the  wave-normal  angle  varies  from  66  deg  at  low  energy  to 
71  deg  at  high  energy  with  a  corresponding  variation  in  ray  angle  from 
156  to  149  deg. 

In  summary,  the  beam  may  send  out  electromagnetic  waves  in  many 
directions  depending  on  beam  energy  and  modulation  frequency;  similarly, 
the  radiation  energy  may  flow  in  many  directions.  Frequencies  below 
100  Hz  and  beam  energies  below  10  keV  are  the  only  regions  in  which 
the  directions  of  the  wave  normal  and  the  ray  are  almost  equal.  Night¬ 
time  emission  at  low  altitude  and  daytime  emission  at  high  altitude 
are  exceptions  as  noted  above.  Rules  for  emission  into  higher  harmon¬ 
ics,  which  are  not  presented  here,  appear  to  be  much  more  complicated. 

EFFECT  OF  VARYING  PITCH  ANGLES 

Figures  39  through  46  show  the  effect  of  varying  the  pitch  angle 
on  both  day  and  night  emission  at  200  km,  A  comparison  of  these  plots 
with  those  for  emission  at  45  deg  (Figs.  9  through  16),  demonstrates 
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that  changing  the  pitch  angle  does  not  greatly  affect  the  peak  values 
of  the  structure  in  the  radiated  power,  but  increasing  the  pitch  angle 
moves  and  stretches  out  the  structures  toward  higher  beam  energies. 

Thus  beams  of  higher  energies  (greater  than  1  kV)  may  not  radiate  at 
zero  pitch  angle,  but  will  radiate  at  pitch  angles  near  90  deg. 

CHARACTERISTICS  OF  CYCLOTRON  MODES 

Figures  47  and  48  show  the  power  radiated  into  the  regular  (s  ■»  1) 
cyclotron  mode  for  the  daytime  ionosphere.  Figures  49  and  50  show  the 
power  into  the  anomalous  (s  ■  -1)  cyclotron  mode.  Figures  51  through 
54  represent  the  corresponding  nighttime  ionosphere.  The  assumed  pitch 
angle  in  all  these  figures  is  45  deg.  The  plots  show  that  very  little 
power  goes  into  root  2,  except  for  the  nighttime  case  above  1  MHz.  In 
that  region,  the  power  levels  for  regular  cyclotron  emission  are  about 
the  same  as  for  Cerenkov  emission. 

Figure  55  shows  contours  of  the  index  of  refraction  at  resonance 
for  regular  cyclotron  emission  using  root  1  during  the  daytime.  Near 
1  MHz,  the  index  of  refraction  varies  from  900  at  10  eV  to  about  3  at 
10  keV.  Near  10  kHz,  it  ranges  from  200  at  the  low  energy  end  to 
about  20  at  10  keV.  Inside  the  rectangular  region  thus  defined,  the 
power  levels  for  regular  cyclotron  emission  are  about  2  orders  of 
magnitude  less  than  for  Cerenkov  emission. 

For  anomalous  cyclotron  emission  into  root  1,  the  power  levels  on 
the  central  plateau  are  also  about  2  orders  of  magnitude  less  than  for 
Cerenkov  emission.  The  plateau  is  stretched  out  to  higher  energies  at 
reduced  power.  Values  for  the  index  of  refraction  are  comparable  to 
the  case  of  regular  cyclotron  emission. 

The  power  levels  at 'frequencies  less  than  60  Hz  are  much  larger 
than  those  for  Cerenkov  emission.  However,  the  resonant  index  of  re¬ 
fraction  is  also  much  larger — about  10  .  As  the  wave  travels  from 

the  ionosphere  to  the  ground  it  enters  a  region  where  the  index  of 
refraction  is  one,  and  therefore  it  has  a  high  probability  of  total 
internal  reflection.  Apparently,  the  cyclotron  modes  will  be  important 
for  ionosphere  to  ground  communication  only  for  frequencies  above  sev¬ 
eral  hundred  kilohertz  at  moderate  electron  beam  energies,  or  above 
10  kHz  for  near-relativistic  electron  beams. 


at  200  km  (root  1,  future  gun 


anomalous  cyclotron  mode  for  daytime  ionosphere 
200  km  (root  1,  future  gun). 


o'oooof; 


noun  i 


to  regular  cyclotron  mode  for 
at  200  km  (root  1 ,  future  gun 


30000.0 


>to  regular  cyclotron  mode  for  nighttime  ionosphere 
at  200  km  (root  2,  future  gun) . 


ooooot; 


clotron  mode  for  nighttime  ionos 


(Sioor  OOOT  OOT  OT 


ro  too 


oooooc 


Resonant  value  of  index  of  refraction  for  regular  cyclotron 
mnHp  at  ?00  km  during  davtiroe  (root  1,  future  gun). 


EFFECT  OF  RETURN  CURRENT 


Presently,  no  adequate  theory  predicts  the  distribution  of  current 
that  returns  to  a  rocket  or  spacecraft  to  neutralize  the  positive  charge 
built  up  as  an  electron  beam  leaves  an  electron  gun,  Harker  has  sug¬ 
gested  that  a  worst-case  estimate  of  the  effect  on  the  radiation  could 
be  modeled  within  the  framework  of  his  theory  by  supposing  that  the 
parallel  components  of  the  emitted  current  are  completely  canceled  by 
the  return  current.  The  expression  for  the  radiated  power  [Eqs.  (3) 
and  (4) ]  is  then  evaluated  using  only  the  terms  that  involve  the 
perpendicular  component  of  the  velocity.  As  might  be  imagined,  the 
radiated  power  is  drastically  reduced  for  pitch  angles  close  to  0  deg. 

At  45  deg  (as  Figs.  56  and  57  show),  for  emission  into  the  nighttime 
ionosphere  at  200  km,  the  root  1  levels  are  reduced  by  9  orders  of 
magnitude  on  the  central  plateau  but  by  less  than  a  factor  of  2  below 
60  Hz.  These  figures  should  be  compared  with  those  for  no  return 
current — Figs.  14  and  16.  The  power  levels  on  root  2  are  reduced  by 
1  order  of  magnitude  in  the  central  region  and  by  5  or  6  orders  of 
magnitude  above  1  MHz.  Daytime  results  are  similar. 

£ 

Although  this  is  a  worst-case  estimate  in  the  context  of  Harker 's 
theory,  it  is  not  based  on  a  self-consistent  calculation  of  the  elec¬ 
tric  fields  built  up  around  a  spacecraft  or  rocket  or  within  the  elec¬ 
tron  beam  itself.  To  our  knowledge,  such  a  calculation  has  not  been 
done.  Some  researchers  believe  that  any  return  current  will  involve 
electrons  moving  at  much  lower  velocities  along  various  trajectories 
and  very  little  cancellation  of  the  parallel  current  will  occur.  (See, 
for  example,  Lavergnat  and  Lehner.^) 

The  collisional  slowing  and  broadening  of  the  electron  beam  is 
related  to  the  return  current.  At  100  km,  the  collisional  effects  may 
be  severe,  yielding  a  radiated  power  lower  than  that  given  in  this 
report.  A  detailed  radiation  emission  and  propagation  calculation 
should  include  electron  collisions  for  beam  injection  at  low  altitudes. 


CONSTANT  BEAM  CURRENT 


In  all  the  calculations  presented  here,  the  power  in  the  electron 
beam  was  a  constant  30  kW.  A  subset  of  the  parameter  study  was  also 
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performed  for  constant  beam  current  at  100  A.  Figures  58  through  61 
show  the  results  for  Cerenkov  emission  with  a  pitch  angle  of  45  deg 
into  the  nighttime  ionosphere  at  200  km.  In  Figs.  58  through  61,  the 
beam  power  is  30  kW  at  beam  energy  of  0.3  keV.  The  beam  power  is,  of 
course,  higher  or  lower  than  30  kW  at  higher  or  lower  energies,  re¬ 
spectively.  The  peak  radiated  power  at  higher  energies  is,  therefore, 
higher  than  that  shown  in  Figs.  13  through  16,  and  lower  at  lower 
energies,  for  emission  at  constant  beam  power,  The  constant  current 
calculation  has  an  additional  feature — it  gives  the  sum  of  the  radi¬ 
ated  power  into  the  frequencies  and  first  five  odd  harmonics.  For  that 
reason  the  forbidden  region  for  root  1  below  1  kHz  appears  to  be 
narrower  than  the  calculations  at  constant  beam  power  would  indicate 
and.  an  additional  allowed  band  is  seen  in  root  2  above  100  kHz.  The 
higher  harmonics  in  the  square  wave  electron  beam  current  evidently 
excited  allowed  electromagnetic  modes  although  the  fundamental  was  in 
a  forbidden  region.  (Note  that  the  plots  show  the  beam  modulation 
frequency,  not  the  radiation  frequency,  which  in  the  band  from  100  to 
4  kHz,  is  above  1  MHz.) 
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IV.  BERT- I 


Plans  call  for  the  BERT-I  rocket  to  be  launched  at  night  from 
White  Sands,  New  Mexico  and  to  travel  in  a  northerly  direction  reach¬ 
ing  an  apogee  of  230  km.  An  electron  beam  will  be  ejected  perpendic¬ 
ular  to  the  rocket  axis,  and  the  rocket  will  be  spinning  during  the 
ejection.  If  the  rocket  is  traveling  due  north  during  the  beam  ejec¬ 
tion,  pitch  angles  from  about  25  to  155  deg  will  be  swept  out. 

The  electron  gun  can  eject  currents  of  0.02,  0.2,  2.0,  and  20  mA 

g 

at  energies  of  0.5,  1,  and  2  keV.  It  may  be  possible  to  modulate  the 
beam  using  a  square  wave  up  to  a  frequency  of  1  kHz.  One  kilohertz 
is  below  the  central  plateau  for  root  1  emission,  so  to  estimate  the 
maximum  possible  power  going  into  higher  harmonics,  the  radiation  up 
to  the  tenth  harmonic  is  kept  in  the  BERT-I  calculation.  Table  2 
summarizes  the  ionospheric  parameters  assumed  for  the  purpose  of  the 
calculation. 


TABLE  2. 

Ionospheric  Parameters  for  BERT-I  Calculation. 


Ion 

Electron 

Ion 

Lower 

Electron 

Gyro- 

Plasma 

Plasma 

Hybrid 

Gyrofrequency 

frequency 

Frequency 

Frequency 

Frequency 

(Hz) 

(Hz) 

Time 

(Hz) 

(Hz) 

(Hz) 

1.36  x  106 

46.1 

Night 

1.77  x  io6 

1.05  x  io4 

7.92  x  103 

Plots  were  generated  for  pitch  angles  of  30,  60,  and  89.9  deg. 
The  extremes,  30  and  89.9  deg,  are  shown  in  Figs.  62  through  69.  Two 
regions  are  of  interest:  one  at  frequencies  lower  than  50  Hz,  the 
other  near  1  kHz.  Figures  62  through  69  show  the  power  radiated  into 
root  1  below  50  Hz.  The  values  are  quite  small — between  10  ^  and 
10  ^  W/m.  The  radiation  into  root  2  is  somewhat  higher  in  this 
region,  and  more  power  goes  into  higher  harmonics  than  into  the 
fundamental. 
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Figures  66  through  69  show  power  into  root  1  near  1  kHz.  In  that 
region  no  power  at  all  goes  into  the  fundamental;  the  fifth  harmonic 
is  the  lowest  frequency  radiated.  The  power  into  root  2  is  6  orders 
of  magnitude  smaller  than  that  into  the  fifth  harmonic  of  root  1  and 
because  it  is  so  small  it  is  not  shown. 

The  figures  confirm  that  the  radiated  power  levels  are  small — 
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on  the  order  of  10  to  10  W/m.  Analysis  suggests  that  even  so,  the 
beam  has  high  enough  gain  (in  the  ideal  case  described  by  Harker  and 
Banks  )  that  radiation  might  be  detectable  at  the  ground.  For  that 
reason,  an  experiment  utilizing  the  BERT-I  beam  might  be  successful; 
however,  there  is  a  high  risk  of  failure. 


V.  VCAP/ SPACELAB  2 


Spacelab  2,  earmarked  for  mid-1985,  will  be  conducted  at  an  alti¬ 
tude  of  350  km.  The  inclination  of  its  orbit  will  be  50  deg,  which 
will  take  the  spacecraft  into  high,  middle,  and  low  geomagnetic  lati¬ 
tudes.  Stanford  University  will  use  the  VCAP  electron  gun  aboard  this 
flight  and  will  launch  beams  at  various  pitch  angles  and  modulation 

frequencies.  The  modulation  will  be  square  wave;  the  frequencies  will 
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range  from  a  few  hertz  to  several  megahertz.  The  gun  currents  will 
be  fixed  at  50  and  100  mA.  The  gun  voltages  may  reach  1  kV.  The  duty 
cycle  may  also  be  varied. 

Table  3  summarizes  the  ionospheric  parameters  PSR  assumed  in  this 
study.  The  analysis  was  geared  toward  a  specific  position  at  midlati¬ 
tude,  above  White  Sands.  A  50  percent  duty  cycle  (equal  on  and  off 
time)  was  assumed  for  the  pulses  and  only  radiation  into  the  Cerenkov 
mode  was  calculated.  The  gun  current  was  assumed  to  be  100  mA  and  the 
first  10  harmonics  were  summed.  Pitch  angles  of  0.1,  30,  60,  and 
89.9  deg  were  run.  Figures  70  through  76  and  77  through  84  show  con¬ 
tours  for  daytime  and  nighttime  emission,  respectively. 

The  figures  reveal  power  levels  up  to  0.1  W/m  on  the  central 
plateau  of  root  1,  much  lower  levels  at  ELF,  and  lower  levels  into 
root  2  (except  above  1  MHz  for  nighttime  emission) .  As  in  the  -case  of 
the  BERT-I  calculation,  the  power  goes  into  higher  harmonics  for  modu¬ 
lation  frequencies  between  6  and  60  kHz. 

TABLE  3. 

Ionospheric  Parameters  for  VCAP/Spacelab  2  Calculation. 

Ion  Electron  Ion  Lower 

Electron  Gyro-  Plasma  Plasma  Hybrid 

Gyrofrequency  frequency  Frequency  Frequency  Frequency 

(Hz)  (Hz)  Time  (Hz)  (Hz)  (Hz) 


Root  1 

No  Return  Current 


Return 


radiated  power  (Cerenkov 
angle  60  deg  (root  1 
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VI.  SEPAC 


The  SEPAC  experiment  is  tentatively  scheduled  to  fly  on  shuttle 
missions  in  August  1986,  in  December  1989,  and  December  1991.  The 
SEPAC  includes  an  electron  gun  that  presently  can  be  modulated  at 
frequencies  from  0.1  to  8.9  kHz.  The  gun  is  capable  of  12.5  kW  emis¬ 
sion  at  highest  voltage,  which  is  7.8  kV.  The  gun  perveance  is 
3/2 

2.5  x  10  A/V  ,  which  limits  the  current  possible  for  a  given 
voltage.  An  altitude  of  275  km  and  a  57  deg  orbit  inclination  are 
planned  for  the  August  1986  flight. ^  As  in  the  cases  of  BERT-I  and 
Spacelab  2,  we  have  assumed  ionospheric  parameters  appropriate  to  a 
White  Sands  overflight. 

Table  4  summarizes  the  calculation  parameters.  In  this  set  of 
runs,  20  harmonics  were  kept,  the  gun  perveance  determined  the  cur¬ 
rent,  and  5  pitch  angles  were  considered.  Three  pitch  angles  are 
presented  here:  0.1,  45,  and  89.9  deg.  Figures  85  through  89  show 
the  power  radiated  into  the  Cerenkov  mode  for  daytime  emission,  and 
Eigs.  90  through  94  show  nighttime  emission.  The  perveance  reduces 
beam  power  and  current  at  low  beam  energies,  therefore,  the  radiated 
power  decreases  as  the  beam  energy  decreases.  This  is  opposite  to 
the  general  tendencies  reported  for  BERT-I  and  Spacelab  2. 

When  the  beam  energy  is  greater  than  1  kV,  the  power  into  root  1 
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is  about  0.1  W/m  and  that  into  root  2  is  about  10  W/m.  For  modu¬ 
lation  frequencies  above  1/2  kHz,  much  of  the  power  goes  into  higher 
harmonics. 

TABLE  4. 

Ionospheric  Parameters  for  SEPAC  Calculation. 


Ion  Electron  Ion  Lower 

Electron  Gyro-  Plasma  Plasma  Hybrid 

Gyrofrequency  frequency  Frequency  Frequency  Frequency 

(Hz)  (Hz)  Time  (Hz)  (Hz)  (Hz) 


1.36  x  IQ6 
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Night  3.36  x  io6  1.96  x  iq4 
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